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Abstract 



■ We propose a model for the quark masses and mixings based on an ^4 family 



symmetry. Three scalar SU (2) doublets form a triplet of A4. The three left-handed- 
I quark SU{2) doublets are also united in a triplet of A4. The right-handed quarks 

are singlets of A4. The A4-symmetric scalar potential leads to a vacuum in which 
two of the three scalar SU (2) doublets have expectation values with equal moduli. 
Our model makes an excellent fit of the observed iKib/^fel- The symmetry CP is 
' respected in the charged gauge interactions of the quarks. 

o ■ 

> ■ 

1 Introduction 

H , 

5r ! In the standard SU{2) x U{1) gauge model for the electroweak interactions, the Yukawa 
couphngs of the fermions to the unique scalar SU{2) doublet of the model are completely 
arbitrary — as a matter of fact, those couplings make up almost all the free parameters of 
the Standard Model (SM). As a consequence of this fact, the SM leaves the quark masses 
and mixings unpredicted; although the model accounts for the existence of quark masses 
and mixings, their actual values remain arbitrary in the context of the SM. 

Several more complex models have tried to overcome this shortcoming of the SM. 
However, most of those models are in reality Ansdtze: instead of deriving the structure of 
the Yukawa couplings from some underlying symmetry of a self-consistent gauge theory, 
they simply assume the Yukawa couplings to have some aesthetically appealing pattern 
or texture. A model should instead rely on some flavour (family) symmetry. 
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If the flavour symmetry is Abelian, then all its irreducible representations (irreps) are 
one-dimensional and the symmetry can at most force some Yukawa couplings to vanish^ 
A non-Abelian flavour symmetry can also force non-vanishing Yukawa couplings to be 
interrelated among themselves through definite Clebsch-Gordan factors. Since there are 
three families of quarks, a most desirable non-Abelian flavour symmetry ought to have 
three-dimensional irreps, in order to achieve full unification of the three generations and, 
hence, to achieve a minimal number of independent Yukawa couplings. The smallest 
discrete group with a three-dimensional irrep is the group of the even permutations 
of four objects. The group A4 has 12 group elements, one triplet irrep 3 and three 
inequivalent singlet irreps 1, 1', and 1" (the 1 is the trivial representation, the irreps 1' 
and 1" are complex-conjugate of each other). This group has, in the last few years, been 
used in many models for the lepton masses and mixings [3J. It has also been used in 
models for the quark sector [4j, or for both quarks and leptons simultaneously [5]. 

In this paper we suggest a model for the Yukawa couplings of the quarks based on 
an ^4 family symmetry. Our model has three scalar gauge-S't/(2) doublets united in a 
3 of A4. The left-handed-quark SU{2) doublets are also united in a 3 of A4. In each 
electric-charge sector, the three right-handed quarks are in a 1 © 1' © 1" of A4. Thus, 
our model achieves a high degree of simplicity and, even, uniqueness, because it treats 
the three families of quarks in the same way and it treats both electric-charge sectors in 
the same way. Furthermore, our model does not require ^4 to be broken anywhere in the 
Lagrangian, not even through soft terms — the breaking of A4 is solely spontaneous. This, 
too, adds to the simplicity of the model. 

Surprisingly, our model is able to predict the mixing parameter iKife/Kbl (where V is 
the quark mixing, or CKM, matrix) fully right: it predicts |Vufe/K:f,| ~ 0.088, in agreement 
with the usual averages of the various phenomenological analyses. On the other hand, 
our model also leads to a null violation of the discrete symmetry CP in the charged gauge 
interactions of the quarks; thus, the observed CP violation, for instance in K^-K^ mixing, 
or in decays, should in the context of our model be explained through scalar-mediated 
interactions, including flavour-changing neutral Yukawa interactions. 

The plan of our paper is the following. In section [2] we derive the form of the quark 
Yukawa- coupling matrices. In section [3] we study the scalar potential and the ensuing 
vacuum. In section H] we write down the quark mass matrices and demonstrate that in 
our model there is no CP violation in the CKM matrix. In section [5] we explain the 
method that we used in the numerical analysis and give some fits and results. A short 
summary is provided in section [61 

2 The Yukawa couplings 

The gauge symmetry of the model is SU{2) x U{1). There are three scalar SU (2) doublets 
(pj [j = 1, 2, 3) with hypercharge 1/2. They form a triplet 3 of the flavour symmetry ^4. 
There are three left-handed-quark SU{2) doublets Q^j with hypercharge 1/6. They are 

^The converse of this statement also holds: any pattern of vanishing Yukawa couplings may be enforced 
by an Abelian flavour symmetry with an adequate spectrum of scalars [I] . 

useful list of all the discrete non-Abelian groups with 31 or less group elements is provided in [2]. 
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united in another 3 of A4. There are three right-handed-quark SU{2) singlets n^j with 
hypercharge —1/3 and three right-handed- quark SU{2) singlets pjij with hypercharge 2/3. 
The uri and pri are y44-invariant (they are I's of A^), the nji2 and pr2 are l"s of A4, 
and the hr^ and pr^ are l'"s of A4. This means that there exist two non-commuting 
transformations Ti and T2, 




Qli, 
* uj'^riRs, 



2, -V3, 

Ql2, Ql3 —Ql3, 



UJ 



(2) 



exp(2i7r/3) = 



under which the Lagrangian is invariant. In equation ( 
(-1 +^73) /2. 

The 3 is a real representation of ^4. Indeed, the scalar SU{2) doublets with hyper- 
charge — 1/2 

0j = iT2(p* (3) 

transform under Ti and T2 in exactly the same way as the 0j, as is obvious from equa- 
tions (H]) and (E]). 

Given this spectrum of fields and their transformation laws under both the gauge 
symmetry SU{2) x f/(l) and the flavour symmetry ^4, the quark Yukawa Lagrangian is 



-^^Yukawa = -Vl (Q Ll4>l + <5l202 + QlS^s) ^ri 



-y2 [Qli4>1 + UjQl24>2 + QL3h) nR2 
-y3 (Qli4>1 + Uj'^Ql2<P2 + t^QL303 ) ^R^i 



-y4 {Qli4>1 + <5l202 + Ql34>-3) PRI 

-1/5 (Ql101 + UjQL2(p2 + ^^^Ql303) PR2 

-2/6 (Ql101 + C^^Ql202 + ^^Qish) PR3 + H.C., 



the six Yukawa couplings ?/i-6 being in general complex. 
The scalar doublets 



(4) 



(5) 



are assumed to have vacuum expectation values (VEVs) 



















V2e''^\ 











^3, 



(6) 



where ^1,2,3 are, without loss of generality, real and non-negative. Since Qlj = (plj, nij) 
the quark mass matrices, defined through 



>^^mass = - ( ^Ll '^L2 '^L3 ) M„ 




PLI PL2 PL3 ) M, 



( Pri 

PR2 
V PR3 



H.C., (7) 



are 



Mr, 



D 



( yivi y2Vi 
yiV2 ujy2V2 

V ym ^^y2V3 
yiVi y^vi 



D* 



y4V2 
y4V3 



^ybV2 

^^y5V3 



ysvi 

^'^y2,V2 

^ysvs 

y&vi 
^'^y&v2 



where 



D = diaer f e 



-ia/2 iPl2 



(9) 



(10) 



The quark mass matrices in equations ([8]) and ([9]) are identical to the mass matrix for 
the charged leptons in the original model of Ma and Rajasekaran (MR) [6]. Indeed, 
we have adopted in our model the same A/^ representations for the quarks as MR did in 
their paper for the leptons. The crucial difference between the two models is that, while 
MR have assumed a vacuum state characterized by v\ = V2= and a = /3 = 0, we shall 
demonstrate the existence of, and employ, a vacuum state with different features. 

Let the unitary matrices f/^'^ satisfy 



y2Vi 



yivi 

yiV2 ujy2V2 ^^yzV2 

2„ 



y^vi 

,2„ 



u 



Pt 



yiV3 uj^y2V3 
I yiVi y^vi 

yiV2 



^y'iVz 

yevi 
ujy^V2 uj'^y&V2 

2„ 



jjri 



diag {md,ms,'mb) , 
diag {mu,mc,mt) . 



(12) 



V y^vs uj^y^v^ ujyQV^ 
Then, the quark mixing (CKM) matrix is 

V 



Ul^D'Ul. 



(13) 



One may absorb the phases of 1/1,2,3 in the overall phases of the three rows of U^, 
and similarly absorb the phases of 1/4,5,6 in the matrix U^. Those six phases are therefore 
unphysical. Thus, this model for the quark masses and mixings has ten parameters: the 
two phases a and (3 in the diagonal matrix D^, and the eight real quantities 



(14) 



As we shall see in the next section, the 744-symmetric scalar potential is so constrained 
that these ten parameters reduce to only eight. 





m 




^3 




m 




ye 


V2 


^1 


Il/l|^3, 








, \y4^3, 














yi 




yi 




1/4 




2/4 


^^3 


V3 



3 The scalar potential 

The most general renormalizable scalar potential invariant under the symmetry A4 is 

V = [(f)\<Pl + 0202 + 0^03) + Al (0l0l + 0^02 + 0^03 
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+ A2 



(402) + (402) (4^3) + {4h) (40l) 



+ (A3 - A2) 
A4 



+ 



0l02) {4<t>l) + (0^3) (402) + (0i0l) (0l0^ 
(0l02)'+(0^03)'+(0l0] 



+ H.C. 



(15) 



where /i and A1-4 are real. The phase e is arbitrary. 
We define v = 



f 1 + f 2 + and 
O3 



e-/3, 
e — a, 
e + a + /3. 



(16) 
(17) 
(18) 



Then, 



Vo = {0\V\0) 



2| \ 4 . \ f 2 2 . 22, 22 
/if + Aif + A3 I f]^f2 + ^2% + "^3^1 

+ A4 (f if 2 COS 6*3 + f 2"^! COS 6*1 + V^vl COS 6^2 

The equations for vacuum stabihty are 
dVo 



dvi 

da 

d(3 














/X + 2Aif ^ + A3 (f 2 + fs) + A4 (f 2 cos 6^3 + f 3 cos 62 
H + 2Aif ^ + A3 (f 1 + f 3) + A4 [yl cos 6^3 + f 3 cos 



1 ' 



/X + 2Aif ^ + A3 {v\ + f 2) + A4 (f 1 cos 62 + f 2 '^os 6*1) 
A4f 1 (— f 2 sin 6^3 + f 3 sin 6*2) , 
A4f 2 (~^i sin 6*3 + f| sin 6^1 



(19) 

(20) 

(21) 

(22) 

(23) 
(24) 



We reject possible solutions to these equations in which one of the VEVs vanishes, and 
also solutions in which fi = f2 = ^3^1 Then, equations (12^ and (1211) yield 



sin 9^ = kv. 



(25) 



where /c is a real constant with dimension M ^. Subtracting equations fl2Tl) and fl22|) from 
equation (I2U1) . one obtains 



(f I — f 1 ) A3 + [(f I — f 1) cos 6*3 + f I (cos ^^2 — COS ^i)] A4 = 0, 
(f I — vD A3 + [(f I — f 1) COS ^2 + vl (cos 6'3 — COS 61)] A4 = 0. 



(26) 



Equations (126!) constitute a Cramer system for A3 and A4. The Cramer determinant must 
vanish and one hence obtains 



aj cos 6'j = 0, 



(27) 



•^Notice that Ma and Rajasekaran opted precisely for vi = V2 = ^3 in their model p], which is otherwise 
similar to ours. 
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where 



2„,2 



2 2 



01 = "Wg — f 2 + tllf 

_ 4 4,2 2 2 2 

02 = Vi - + V2V3 - V1V2, 

_ 4 4,2 2 2 2 

as = 1-2 - f 1 + Wifg - 1-2^3 



satisfy 



Equation (p7|) together with equation (!25l) imply 



-A aWl - asVl - ^^^^2, aaV ^ ~ ^^^3 



where [7] 



A (a, b, c) = -a^ -h^-c^ + 2a%'' + 2a^^ + 2hW 
has five solutions: 



(28) 
(29) 
(30) 



(31) 



(32) 
(33) 



Equation 

1. A; = 0, 

2. P = A(t;f,t;i,^i)/(4t^^2M), 

3. Vl = V2, 

4. Vl = V3, 

5. f2 = V3. 

It is easy to explore in detail solutions 1 and 2 and to show that they require A3 = ±A4 
(solution 1 furthermore needs e = 0). Thus, those solutions require a non-trivial constraint 
on the parameters of the potential, and that constraint is in general unstable under 
renormalization. Those solutions should therefore be discarded. The remaining solutions 
3-5 are equivalent, since the three scalar doublets 0j form a triplet of ^44. We shall use 
for definiteness solution 3: Vi = V2 and 9i = 62, i.e. a = (3. Equations fl2(l - fl2^ then 
reduce to only three equations. 



= /i + 2Ai (2vl + vl) + As {vl + vl) + A4 
= /i + 2Ai {2vl + vl) + 2Xsvf + 2X4vf cos (e 
= t;^ sin (e + 2a) — t>3 sin (e — a). 



vl cos (e + 2a) + vl cos {e — a) , (34) 



a) 



(35) 
(36) 



which determine the three quantities Vi, v^, and a. 
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4 The mass matrices and CP conservation 



We saw in the previous section that, just as we had advertised at the end of section [21 
consideration of the most general A4-invariant scalar potential actually reduces the ten 
parameters of our model to only eight, since V2 = fi and (3 = a. The quark mass matrices 
of our model are therefore 



M„ 



Mr, 



diag (e-*"/2^ e*"/2^ l) 



diag (e'"/^ e-'"/2, l) 



a b 
a ub 

ra uj'^rb ujrc 

( f 9 h 

f ujg uj'^h 

\ rf uo'^rg ujrh 



(37) 
(38) 



where r = v-^/vi and a, 6, c, /, g, and h are real and positive. Then, 



/ 



y e 

X rye 

ia/2 



ia/2 



X ye rye 



(39) 



(40) 



ye' 

\ ry*e*"/^ rye 
z w* 

where x = a? + b'^ + (? and z = p + g"^ + are real, while y = + ub"^ + tu^c^ and 
w = p + uig"^ + u'^h'^ are complex. Now, computing the commutator of Hp and Hn one 
finds that it is of the form 



2 

r X 

-ia/2 

2 

r z 



[Hp, Hn 



f -n 
n 
\ m* m 




(41) 



hence det [Hp, Hn] = 0. Therefore, in this model there is no CP violation in the quark 
mixing matrix, i.e. the Jarlskog observable J [8] vanishes. 



5 Numerical procedure and results 

We have performed a global analysis of the the quark mass matrices given in the 
previous section — equations (1371) and (!38l) — by employing the downhill simplex method [9] . 
Table [T] specifies in its first two columns the observable quantities in the form 

Oi = Oi±ai, (42) 

where Oj is the experimental mean value of Oj and cxj is the square root of its variance. 
The index i = 1, . . . ,9 labels the nine observables given in tabled! Writing x for the set of 
the eight parameters of our model (a, b, c, f, g, h, r, and a), and Pi (x) for the resulting 
predictions for each of the observables, one constructs the function 

(x) = E 

i=l 



P U) 



o, 



(43) 
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Observable 


Experimental value 


Model prediction 


Pull 


rrid [MeV] 


5±2 


4.977 


-1.2 X 10" 


-2 


rus [MeV] 


95 ±25 


90.545 


-1.8 X 10' 


-1 


rrib [MeV] 


4200 ± 70 


4200.79 


±1.1 X 10' 


-2 


rriu [MeV] 


2.25 ±0.75 


2.250 


±1.9 X 10" 


-5 


rric [MeV] 


1250 ±90 


1250.498 


±5.5 X 10" 


-3 


mt [GeV] 


172.5 ±2.7 


172.497 


-1.2 X 10' 


-3 


sin 


0.2243 ±0.0016 


0.22431 


±8.1 X 10' 


-3 


sin 023 


0.0413 ±0.0015 


0.04139 


±6.0 X 10' 


-2 


sin 6li3 


0.0037 ±0.0005 


0.003627 


-1.5 X 10' 


-1 



Table 1: Experimental data and result of our best fit. The experimental data (average 
values and error bars) used in our numerical analysis are given in the second column. 
The data on the quark masses have been taken from p^. The light-quark masses are 
renormalized at a common scale ~ 2 GeV. The charm- and bottom-quark masses are 
running masses in the MS scheme. The data on the quark mixing angles have been taken 
from pT] . The third column displays the values Pi predicted by our model when the values 
of its parameters are those in equations flUl) . The fourth column shows the number of 
standard deviations from the mean values, (^Pi — Oij /ai, computed using the data from 
the second column. The value = 0.057 is the sum of the squares of the numbers in the 
fourth column and is dominated by the pulls of and sin ^13. 

The global minimum of represents the best possible fit of the model predictions to the 
experimental data. 

We found an excellent fit, with = 0.057, of our model to the nine input data 
specified in table [H The input parameters of the fit are 



a 


= 40.75189 MeV, 


b 


= 87.78761 MeV, 


c 


= 2.347665 MeV, 


f 


= 3941.127 MeV, 


9 


= 515.0460 MeV, 


h 


= 1.060808 MeV, 


r 


= 43.37746, 


a 


= 0.2251660 radians. 



(44) 



Other details of the fit are given in the third and fourth columns of table [H 

Notice that the observables in table [1] do not include the phase 6 of the quark mixing 
matrix. Indeed, since we already know that in our model that matrix is real, we are of 
course unable to fit for any 5 7^ 0, tt. In our model the explanation of all the observed 
CP violation through a sole phase in the quark mixing matrix gets spoiled: other sources 
of CP violation must be found, and fitted to the observed data. In our model the most 
likely source of CP violation will be the interactions mediated by both charged and neutral 
scalars, since ours is a three-Higgs-doublet model featuring, in particular, flavour- changing 
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interactions mediated by neutral scalars. 

In order to test the variation of as a function of the value Oj of an observable 



2 



quantity Oj, we substitute in the expression for (x) the term Pj (x) — Oj / [oif by 



a term 



Pi fx) - O, 



/ (0.01 Ojj . The small error assigned to Oj in the denominator of 

this term guarantees that Oj gets pinned down to the value Oj. 

Figure [1] depicts as a function of r, i.e. of the ratio of VEVs v^/vi. We read off 
from that figure that only for 40 ^ r ^ 52 can good fits be obtained; thus, the range of 
the ratio of VEVs is severely constrained. 

In figure [2] (left panel), the change of under variations of the quark-mixing observ- 
able iKife/Vcbl is shown. There is a pronounced minimum of for 0.08^ iKife/Kbl ^0.09; 
this is in excellent agreement with the value obtained for that observable by the phe- 
nomenological analyses. This remarkable result of our model provides a clear-cut predic- 
tion for |\4f,/\46|. Figure [2] (right panel) gives x^ as a function of |\4d/Vts|. We find in this 
case excellent fits whenever 0.14 < |^td/^s| ~ 0.15. Clearly, this result is correlated to 
our model's prediction for iKife/Kbl? since in our model CP is conserved in quark mixing 
and therefore the CKM matrix is determined by only three parameters. 



6 Conclusions 

In this paper we have proposed a self-consistent model for the quark masses and mixings 
based on a family symmetry ^4. The Yukawa-coupling matrices of our model contain, 
at face value, ten parameters, but, when one considers the y44-symmetric scalar potential 
in detail, one sees that two of those parameters actually disappear. In our model the 
family symmetry ^4 is not broken anywhere in the Lagrangian — its breaking is fully 
spontaneous. The model gives a perfect fit of the observed quark masses and mixing 
parameters, except for the fact that there is no CP violation at all in the CKM matrix. 
The observed CP violation should result in our model from scalar- mediated interactions, 
in particular flavour-changing neutral Yukawa interactions at tree level and also charged- 
scalar-mediated box diagrams at the one-loop level; a detailed study of those interactions 
should be the subject of a separate publication. 

In previous ^4 models by other authors a vacuum characterized by equal VEVs vi = 
V2 = fs has often been employed. In our ^4 model for the quark sector, on the contrary, 
we use f 1 = f 2 while ^3 is some 43 times larger. It is not clear to us whether this vacuum 
may or may not produce a viable model for the lepton sector as well; this will be the 
subject of future investigation. 

Acknowledgement: The work of L.L. was supported by the Portuguese Fundagdo para 
a Ciencia e a Tecnologia through the project U777-Plurianual. 
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